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Abstract. The current evidence about the age of stellar populations in 
elliptical galaxies is reviewed. The case for the bulk of stars in galactic 
spheroids (ellipticals and bulges) having formed at high redshift (z>3) 
is now compelling, both for cluster and well as field galaxies. Whether 
the assembly of ellipticals is deferred to lower redshifts compared to the 
formation of their stars remains controversial, and we mention ongoing 
observational programs that are designed to solve the issue. We present 
preliminary results of a pilot project aimed at ascertain the nature of 
ON ■ extremely red galaxies (ERGs) with near-IR spectroscopy at the VLT. 

S3 

Q ■ 1. The Age of Cluster Ellipticals 

Through the '80s it became normal to start a talk about ellipticals by saying that 

'the classical view of elliptical (and spheroid) formation was that they formed 

stars efficiently in the early universe, and then their stellar population evolved 

k>( | passively after this initial burst'. Then the talk would have continued by saying 

5_i ■ '... however, in recent years there has been growing evidence that this is not the 

true story, ellipticals are rather intermediate age objects, and star formation in 
them has continued for a major fraction of the Hubble time and in some cases 
is still going on today ...'. Actually, such evidence hardly grew, in spite of a 
widespread expectation given the success of hierarchical theories of structure 
formation coupled to the aesthetic appeal of Q = 1 inflationary cosmology. 

^From the turn of the decade through all the '90s the has been instead 
growing empirical evidence that ellipticals are dominated by very old stellar pop- 
ulations. The first breakthrough came from noting the tightness of the color-o" 
relation of ellipticals in the Virgo and Coma clusters (Bower, Lucey, & Ellis 
1992). This demands a high degree of synchronization in the star formation 
history of these ellipticals, that is most naturally accounted for by pushing back 
to early times most of the star formation. Making minimal use of stellar popu- 
lation models, this approach provided for the first time a robust demonstration 
that at least cluster ellipticals are made of very old stars, with the bulk of them 
having formed at z>2. 



The main lines of the Bower et al. argument are as follows. The observed 
color scatter of cluster ellipticals is related to the age dispersion among them by 
the relation: 

8(U -V) = 9{U m V) (t H - t F ) (1) 

where tn and tF are the age of the "oldest" and "youngest" galaxies, respec- 
tively. Here by age one intends the luminosity-weighted age of the stellar popu- 
lations that constitute such galaxies. The time derivative of the color is obtained 
from evolutionary population synthesis models, which give d(U — V)/dt ~ 0.02 
mag/Gyr for t ~ 10. The observed color scatter is 5{U — V) ci 0.04 mag, consis- 
tent with pure observational errors. Hence, one gets tu — ^F ^0.04/0. 02 = 2 Gyr, 
and if the oldest galaxies are 15 Gyr old, the youngest ones ought to be older 
than 13 Gyr, from which Bower et al. conclude they had to form at z>2. If the 
oldest galaxies were instead as young as, say 5 Gyr, then the youngest should be 
older than at least 3 Gyr, which would require a high degree of synchronization 
in their formation, which seems unlikely. 

Evidence in support of the Bower et al. conclusion has greatly expanded, 
to finally become compelling. It is worth mentioning some of the steps of this 
development: the tightness of the fundamental plane relation for ellipticals in 
local clusters (Renzini & Ciotti 1993), the tightness of the color-magnitude re- 
lation for ellipticals in clusters up to z ~ 1 (e.g., Aragon-Salamanca et al. 1993; 
Taylor et al. 1998; Kodama et al. 1998; Stanford, Eisenhardt, & Dickinson 
1998), and the modest shift with increasing redshift in the zero-point of the 
fundamental plane, Mg— a, and color-magnitude relations of cluster ellipticals 
(e.g., Bender et al. 1997; Dickinson 1995, 1997; Ellis et al. 1997; van Dokkum 
et al. 1998; Pahre, Djorgovski, & de Carvalho 1997; Stanford, et al. 1998; 
Kodama et al. 1998). All these studies agree in concluding that most stars in 
ellipticals formed at z^3, though the precise lower limit of the redshift depends 
on the adopted cosmology; it would be more like z>2 in the current standard 
cosmology: {U, A) = (0.3,0.7). 

It is worth emphasizing that all these studies have in common the same 
methodology pioneered by Bower et al. (1992). They focus indeed on the tight- 
ness of some correlation among the global properties of cluster ellipticals, which 
sets a robust constraint on their age dispersion as opposed to an attempt to date 
individual galaxies. Moreover, the move to high redshift offers two fundamental 
advantages. The first advantage is that looking at high z provides the best pos- 
sible way (we should say the way) of removing the age-metallicity degeneracy. 
If spheroids are made of intermediate-age, metal rich stars, they should become 
rapidly bluer and then disappear already at moderate redshift, which is not the 
case (e.g. Kodama &: Arimoto 1997). 

It has also been demonstrated that the scatter about the average Mg— a 
relation (among local z ~ cluster ellipticals) is larger than expected from 
pure observational errors (Colless et al. 1999). Colless et al. interpret this as 
evidence for a substantial age dispersion, but it remain to be seen whether this 
interpretation is consistent with the color scatter remaining constant all the way 
to z ~ 1, as found by Stanford et al. (1998). The observational opportunity of 
studying galaxies at large lookback times makes quite obsolete those attempts 
to find combinations of spectral indices that are aimed at distinguishing between 
age and metallicity effects in nearby galaxies. The obvious drawback of these 



indices (e.g. Hg) is that they are very sensitive to even minor recent episodes of 
star formation, and therefore they do not help to determine the formation epoch 
for the bulk of stars in ellipticals. 

The second advantage of high— z studies is that at high redshift one gains 
more leverage: for given dispersion in some observable one can set tighter and 
tighter limits to the age dispersion. This comes from the color time derivatives 
being larger the younger the population. For example, the derivative d(U— V)/dt 
is ~ 7 times larger at t = 2.5 Gyr than it is at t = 12.5 Gyr (e.g. Maraston 1998), 
and therefore a given dispersion in this rest-frame color translates into a ~ 7 
times tighter constraint on age and therefore on formation redshift. The case is 
effectively illustrated by isolated high redshift ellipticals: Spinrad et al. (1997) 
found a fossil (i.e. passively evolving) elliptical at z = 1.55 for which they infer 
an age of at least 3.5 Gyr, hence a formation redshift in excess of ~ 5. A similar 
case is anticipated by Dunlop (1998), with a galaxy at z = 1.43 and an age of 
3-4 Gyr. At an even higher formation redshift may lie the extremely red galaxy 
in the NICMOS field of the HDF-South, which spectral energy distribution is 
best accounted for by an old, passively evolving population at z — 2 (Stiavelli 
et al. 1999). 

2. Cluster vs Field Ellipticals 

Most of the evidence mentioned in the previous Section refers to cluster ellipti- 
cals. In hierarchical models, clusters form out of the highest peaks in the pri- 
mordial density fluctuations, and cluster ellipticals completing most of their star 
formation at high redshifts could be accommodated in the model (e.g. Kauff- 
mann 1996; Kauffmann & Chariot 1998a). However, these models predict that 
in lower density, field environments, both star formation and merging should be 
appreciably delayed to later times, which offers the opportunity for an observa- 
tional test of the hierarchical merger paradigm. 

The notion of field ellipticals being a less homogeneous family compared to 
their cluster counterparts has been widely entertained, though - once more - the 
direct evidence has been only rarely discussed. Visvanathan & Sandage (1977) 
found cluster and field ellipticals to follow the same color-magnitude relation, 
but Larson, Tinsley, & Caldwell (1980) - using the same database - concluded 
that the scatter about the mean relation is larger in the field than in clusters. 
More recently, a larger scatter in field versus cluster ellipticals was also found for 
the fundamental plane relations by de Carvalho & Djorgovski (1992). However, 
when using absolute luminosities and effective radii at least part of the larger 
scatter among field ellipticals certainly comes from their distances being more 
uncertain than for cluster galaxies. 

Taking advantage of a large sample (~ 1000) of early-type galaxies with 
homogeneously determined Mg2 index and central velocity dispersion, Bernardi 
et al. (1998) have recently compared the Mg2 — a relations (which are distance 
independent!) of cluster and field galaxies. The result is that field and cluster 
ellipticals all follow basically the same Mg2 — a relation. The zero-point offset 
between cluster and field galaxies is 0.007± 0.002 mag, with field galaxies having 
lower values of Mg 2 , a statistically significant, yet very small difference. This is in 



excellent agreement with the offset of 0.009 ±0.002 mag, obtained by Jorgensen 
(1997) using 100 field and 143 cluster galaxies. 

Using the time derivative of the Mg2 index from synthetic stellar popu- 
lations, Bernardi et al. conclude that the age difference between the stellar 
populations of cluster and field early-type galaxies is at most ~ 1 Gyr. The ac- 
tual difference in the mass-weighted age (as opposed to the luminosity-weighted 
age) could be significantly smaller that this. It suffices that a few galaxies have 
undergone a minor star formation event some Gyr ago, with this having taken 
place preferentially among field galaxies. 

3. Spiral Galaxy Bulges 

According to general wisdom, the bulges of spiral galaxies are hard to distinguish 
from elliptical galaxies of the same luminosity, once they are stripped of the 
disk which rotates around them. This is illustrated for example by most bulges 
following the same Mg2 — a and Mg2— luminosity relation as common ellipticals 
(Jablonka, Martin, &: Arimoto 1996). In their sample just ~ 20% of the studied 
bulges appear to be appreciably bluer that the standard relation for ellipticals, 
a sign that recent star formation is required in only a minority of the bulges. 

Closer to us, contrary to several early claims, no evidence exists for an 
intermediate age population for the bulge of M31 (e.g. Renzini 1998; 1999; 
Jablonka et al. 1999). Even more stringently, HST and ground based color- 
magnitude diagrams of Galactic bulge globular clusters and fields indicate a 
uniform old age for our own bulge (Ortolani et al. 1995). Given that we leave in 
a galaxy that is member of a rather loose group, certainly not an exceptionally 
high peak in the primordial density fluctuations, it seems reasonable to conclude 
that star formation in most galactic spheroids, ellipticals and bulges alike, was 
essentially complete at high redshift, no matter whether such spheroids resides 
today in high- or low-density environments. 

4. Do Ellipticals Disappear by z ~ 1? 

In apparent conflict with this conclusion are several claims according to which 
the comoving number density of field elliptical galaxies is rapidly decreasing 
with redshift by z ~ 1 (e.g. Kauffmann, Chariot, & White 1996; Zepf 1997; 
Meneanteau et al. 1998; Franceschini et al. 1998). This is a controversial result. 
Other groups claim such dismissal of high-redshift ellipticals to be premature, 
either because the adopted models for putative high-z ellipticals would be too 
naive (Jimenez et al. 1999), or because they find evidence for the comoving 
number density of ellipticals being constant to at least z ~ 1 (e.g. Totani & 
Yoshii 1998; Benitez et al. 1998; Shade et al. 1999), and possibly as much as to 
z ~ 2 (Broadhurst & Bouwens 1999). 

Part of this discrepancy is likely to arise from the samples so far analyzed 
being rather small, both in terms of the total number of galaxies involved, and 
especially in terms of the (small) number of large scale structures that happen 
to be included in each sample (see below). Part of the discrepancy may also 
arise from the use of different selection criteria for ellipticals. Clearly, these 
differences need to be understood before claiming evidence in support of a late 



assembly of ellipticals (as favored by most current realization of the CDM model 
of hierarchical formation of cosmic structures, e.g. Kauffmann & Chariot 1998a), 
or in support of an early build-up of massive galaxies, resembling the "old- 
fashioned" monolithic collapse models. 

5. Searching for Fossil Ellipticals at High z: Early VLT Observations 
of Extremely Red Galaxies 

The existence and the abundance of fossil ellipticals at z > 1 can be con- 
strained with the selection and study of faint galaxies with the colors expected 
for passively evolving old stellar populations. For instance, a color threshold of 
R— K>5.3 is expected at z^l in the case of a passively evolving, solar metallicity 
stellar population with Z{ OTm > 3, SFR oc exp(—t/r) and r=0.1 Gyr (adopting 
the Bruzual & Chariot 1998 models, and assuming Hq = 50 kms -1 Mpc" 1 and 
Slo = !)• The observed colors of the few known isolated z > 1 ellipticals (see 
Section 1) are indeed consistent with those expected from the spectral synthesis 
models. 

Do such color-candidate ellipticals at z > 1 exist? Observationally, the 
answer is yes: a population of extremely red galaxies (hereafter ERGs) was dis- 
covered with the combination of optical and near-IR imaging (e.g. Elston et al. 
1988; McCarthy et al. 1992; Hu & Ridgway 1994). ERGs are ubiquitously found 
in empty sky fields (Cohen et al. 1999; Thompson et al. 1999), in the vicinity of 
high-z AGN (McCarthy et al. 1992; Hu & Ridgway 1994) and as counterparts of 
faint X-ray (Newsam et al. 1997) and weak radio sources (Spinrad et al. 1997). 
Their typical magnitudes K ~ 18 — 20, and their surface density for R — K > 6 
and K < 19 is 0.039±0.016 ERGs/sq. arcmin (Thompson et al. 1999). Their 
nature is still poorly known, as their faintness makes optical and near-IR spec- 
troscopic observations quite difficult, but they are likely to represent a mixed 
bag class of objects. Indeed, their red colors are consistent with them being pas- 
sively evolving ellipticals, but also with strong dust reddening in a star- forming 
or active galaxy. A dramatic example of this ambiguity is provided by HR10: its 
spectral energy distribution (SED) is consistent with that of an old elliptical at 
z ~ 2.4 (Hu & Ridgway 1994), but near-IR and submm observations proved that 
HR10 is a dusty starburst galaxy at z = 1.44 (Graham & Dey 1996; Cimatti et 
al. 1998; Dey et al. 1999). On the other hand, recent observations suggest that 
the ERG CL 0939+4713B (R - K ~ 7) is an old galaxy at z ~ 1.6, based on 
the detection of the redshifted 4000 A break (Soifer et al. 1999). 

In order to investigate the nature of ERGs and their role in the framework 
of elliptical galaxy formation and evolution, we recently started a project based 
on optical and near-IR imaging and spectroscopy with the ESO NTT 3.5m and 
the VLT 8m telescopes. Near-infrared (JHK) spectroscopy of an incomplete 
"pilot" sample of ERGs was made in April 1999 with the ESO VLT-UT1 (Antu) 
equipped with the IR imager-spectrograph ISAAC (Moorwood et al. 1999). The 
typical integration times were 1-2 hours per spectrum and the resolution was 
~ 500. Nine ERGs with R — K > 5 or I — K > 4.5 were observed with the 
main purpose of finding cases of bona fide ellipticals at z > 1 or, in case of dusty 
starburst/active galaxies, to search for redshifted emission lines (such as Ha at 
z > 0.7) (Cimatti et al. 1999). 



10" 



<< 






tf^|''~" J101950-2220.9 z-1.56 SSPJ2.7 Gyr E(B-V)=0 

I T 




0.5 



1.5 

\(jixa) 



Figure 1. The optical to near-IR SEDs of five ERGs from broad-band photometry 

(filled triangles) and fSAAC near-IR spectroscopy (vertical bars). The spectra are binned 
with 200 A wide bins in order to reach a typical O/N ~10 per bin. The first three SEDs 
from the top are those consistent with both passively evolving ellipticals and no dust ex- 
tinction. Instead, the other two SEDs require substantial dust extinction to be reproduced. 
The continuum and the dotted lines are the best fit spectra obtained respectively with 
dustless, old stellar populations and with dust-obscured constant star formation. The best 
fit parameters (Z, ages and Eg—y) are shown relatively to the best fit of the two cases 
(i.e. dusty or dustless). From bottom to top, the fluxes of the five SEDs are multiplied by 
10 ' ' ' ' respectively. 



Continuum emission was detected in all the targets, but neither strong 
emission lines nor evident continuum breaks were detected in the ISAAC spectra. 
The depth of our spectra (typically Fn m < 2xlCP 16 erg s _1 cm -2 ) would allow 
to detect emission lines such as Ha if as strong as in HR10 (Dey et al. 1999). 
These limits on the Ha emission imply upper limits on the star formation rates 
of individual galaxies that - in absence of dust extinction - range from ~ 7 to ~ 
30/i^q M Q yr _1 (f^o = 1.0). These limits correspond to the SFRs as high as those 
of nearby gas rich spiral galaxies (Kennicutt 1998). However, in case of strong 
dust extinction, the SFRs would increase significantly. For instance, adopting 
the Calzetti (1997) extinction curve, for Eb-v =0.2-0.9, the corresponding 
correction factors of the SFR (as derived from Ha) are ~2-30x. 

In absence of distinctive spectral features, the ISAAC spectrophotometry 
combined with the available optical and near-IR broad-band photometry was 
used to estimate the "spectrophotometric" redshifts of the observed ERGs and 
to study their SEDs. We found that a fraction of the program ERGs have 
SEDs consistent with being passively evolving ellipticals at z ~ 1.5 — 1.8 with 
no dust extinction and rather old ages (~ 1 — 3 Gyr) (see Fig. 1). Their 
rest-frame if-band absolute magnitudes [Mk ~-24.8-t- -25.2) imply luminosi- 
ties L<L* (adopting M^= -25.16 for the local luminosity function of elliptical 
galaxies; Marzke et al. 1998). In addition, deep HST optical imaging shows that 
their morphologies are compact and regular, consistent with such ERGs being 
high-z ellipticals. However, in some cases the SEDs seem to require strong dust 
extinction to be reproduced (see Fig. 1), and in a fraction of these cases the 
morphologies are clearly disturbed (i.e. incompatible with being ellipticals). In 
other cases the required dust may be an artifact of the fitting procedure, in 
which the metallicity of the population is kept fixed to solar. 

On the one hand, our results show the existence of a population of faint 
galaxies that are selectable through their very red colors and which have prop- 
erties consistent with the strict definition of being dustless, old and passively 
evolved spheroidals at z > 1 (see also Spinrad et al. 1997; Soifer et al. 1999). 
Not surprizingly, our observations also indicate that a non-negligible fraction of 
the ERG population consists of galaxies whose SEDs and morphologies are more 
consistent with being high-;?, star-forming dusty systems. 

The sample for this pilot study is clearly too small to draw more quantitative 
conclusions at this stage. We now plan somewhat deeper spectroscopic obser- 
vations of a complete sample, with the aim to assess what fraction of the ERG 
population consists of passively evolving high-z ellipticals and thus to compare 
their observed abundance with that expected from pure luminosity evolution 
models. 



6. Do Bright Galaxies Disappear by z ~ 1? 

Generally, ellipticals have been selected either according to color or morphologi- 
cal selection criteria, or to some combination thereof. However, as one moves to 
high redshift, minor residual star formation may cause virtually fully assembled 
ellipticals to drop out of samples that are selected using color and/or morpholog- 
ical information. To circumvent this limitation, Broadhurst, Ellis, & Glazebrook 
(1992) proposed to select objects only according to their if -band magnitude and 



to construct their redshift distribution. This provides a more robust measure of 
the evolution (if any) of the comoving number density of massive galaxies, and a 
more fundamental check of the merging paradigm. The advantage of a K-band 
selection is that, as opposed to e.g. l?-band selected samples, it favors massive 
galaxies even if currently experiencing very low levels of star formation. 

Following this approach, Kauffmann & Chariot (1998b) estimate that in 
their pure luminosity evolution (PLE) model ~ 60% of the galaxies in a if < 
20 sample should be at z > 1, while only ~ 10% are expected in a standard 
hierarchical merging model. The largest redshift surveys of K < 20 selected 
samples currently available are the Keck samples of Cowie et al. (1996) covering 
~ 26 square arcmin, and Cohen et al. (1998) covering ~ 14 square arcmin. 
Both spectroscopic samples are rather incomplete, for example, the Cohen et 
al. sample includes 195 objects, among which 24 are stars, 21 are galaxies 
at z > 1, and optical spectroscopy failed to provide a redshift for another 32 
objects. The vast majority of these latter objects are likely to be galaxies at 
z>l, in which major spectral features have moved to the near IR. All in all, in 
this sample ~ 53/171 ~ 30% of galaxies are likely to be at z > 1. This empirical 
value falls (ironically enough) just midway between the prediction of the two 
competing models. Most importantly, Cohen et al. emphasize that ~ 50% 
of their galaxies fall in just five 'redshift peaks' most likely due to clustering. 
Therefore, a much larger sample is required to dispose of a 'fair' representation 
of the galaxy population, one that is not subject to statistical fluctuations in 
the number of included large scale structures. 

The conclusion is that the two Keck redshift survey of K -limited samples 
are so far insufficient to conclude the experiment envisaged by Broadhurst, El- 
lis, & Glazebrook (1992). Big samples are required, distributed over several 
different lines of sight so as to average out the effects of clustering and large 
scale structures. A contribution in this direction will come from an observing 
program with FORS1, FORS2 and ISAAC at the VLT (PI A. Cimatti), aimed 
at gathering redshifts of two K < 20 samples over two separate lines of sight, 
for a total of ~ 500 galaxies. 

7. Speculations 

In the previous sections it has been documented that compelling evidence exists 
for the bulk of stars in galactic spheroids being very old, i.e. formed at redshifts 
beyond ~ 3, and possibly even much beyond this value. This applies to ellipticals 
and bulges alike, in clusters as well as in the lower density regions still inhabited 
by spheroids, including the bulge of our own Galaxy. This is what was expected 
(actually postulated) in the monolithic collapse scenario, while it is at variance 
with most existing realizations of the hierarchical merging scenario. 

The fact that spheroids are made of old stars does not necessarily invalidate 
the hierarchical merging paradigm, which actually offers a still unique descrip- 
tion on how large galaxies could have been assembled. In an effort to comply 
with the observations, hierarchical models should be tuned to mimic the mono- 
lithic model as much as possible, primarily by pushing most of the action back to 
an earlier cosmological epoch. This is certainly favored by low-fi + A cosmolo- 
gies, but will not suffice. As universally recognized, the pitfall of scmianalytical 



realizations of the hierarchical paradigm is represented by how star formation 
and its feedback are parameterized in the models. The canonical assumption 
has been that most of the star formation takes place in spiral disks, and el- 
lipticals are formed by merging spirals. We think that this widely entertained 
paradigm is far from being proven. All we know about ellipticals demands "that 
star formation, metal enrichment, merging and violent relaxation did not take 
place sequentially; they more likely were all concomitant processes, along with 
supernova heating and radiative cooling in a multiphase ISM" (Renzini 1994). 

With most of the merging taking place at high redshifts, among still mostly 
gaseous components, merging itself would promote the widespread starburst 
activity responsible for the prompt buildup of galactic spheroids (Somerville 
& Primack 1998; Renzini 1999). The natural observational counterparts of 
these events may be represented by the Lyman-break galaxies at z>3 (Stei- 
del et al. 1999), where star formation rates could be in extreme cases as high 
as ~ 1000 MQyr -1 (Dickinson 1998). It remains to be explored whether such 
tuning of the star formation algorithms and of the many parameters of the 
semiempirical models could produce model universes fulfilling all other observa- 
tional constraints. 

In this scenario disks are not primordial to spheroids. They rather develop 
only later, being accreted by those spheroids that happen to be in a permissive 
environment. 
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